library(dplyr) library(plyr)
#mercury data frames DFfl25miA <- as.data.frame(fl25miA) DFfl100miA <- as.data.frame(fl100miA) DFfl25miB <- as.data.frame(fl25miB) DFfl100miB <- as.data.frame(fl100miB) DFla25miA <- as.data.frame(la25miAsw) DFla100miA <- as.data.frame(la100miAsw) DFla25miBA <- as.data.frame(la25miBAsw) DFla100miBA <- as.data.frame(la100miBAsw)
#relevant mrip data landOtimA landOtimB landTtimA landTtimB
#isolate species, biomass, landings, per site - FLorida
DFfl25miA$Whole_weight <- as.numeric(DFfl25miA$Whole_weight)
NAs introduced by coercion
DFfl25miA$TL <- as.numeric(DFfl25miA$TL)
str(DFfl25miA)
'data.frame': 308 obs. of 37 variables:
$ Group : int 121 120 122 122 122 122 122 122 122 122 ...
$ Bottle : chr "TB04010" "TB04009" "TB0400757-D" "TB0400721-D" ...
$ Location.x : chr "TBM" "TBX" "TB" "TB" ...
$ Field_no : chr "TBM04110906" "TBX04062101?" "" "" ...
$ Spec_no : int NA NA 1 1 2 3 4 2 1 1 ...
$ FIM_spec_no : chr "1" "1" "." "." ...
$ Other_field_no: chr "" "" "1530420040428037" "1000920040423002" ...
$ Date : int 20041104 20040621 20040428 20040423 20040423 20040423 20040423 20040423 20040423 20040421 ...
$ New_date : chr "11/4/2004" "6/21/2004" "4/28/2004" "4/23/2004" ...
$ Year : int 2004 2004 2004 2004 2004 2004 2004 2004 2004 2004 ...
$ Scfname : chr "Rachycentron canadum" "Carcharhinus leucas" "PARALICHTHYS ALBIGUTTA" "EPINEPHELUS MORIO" ...
$ Species : chr "Cobia " "Bull shark " "Gulf flounder " "Red grouper " ...
$ NODCCODE : num 8.84e+09 8.71e+09 8.86e+09 8.84e+09 8.84e+09 ...
$ Sex : chr "F" "F" "F" "U" ...
$ SL : chr "1030" "." "289" "417" ...
$ FL : chr "1122" "." "350" "193" ...
$ TL : num 1279 2800 350 518 535 ...
$ Whole_weight : num NA NA NA 2 2 2 2 2 NA 5 ...
$ Zone : chr "b" "" " " " " ...
$ Grid : chr "." "." "." "." ...
$ Lat_deg : chr "." "." "." "." ...
$ Lat_min : chr "." "." "." "." ...
$ Site : chr "big bayou" "found at gandy ramp" "ruskin" "clearwater" ...
$ Region.x : chr "MW" "MW" "Tampa Bay" "Tampa Bay" ...
$ Hg : num 0.8 1.9 0.16 0.24 0.33 0.28 0.25 0.32 0.6 0.72 ...
$ Lat_Deg : int 27 27 27 27 27 27 27 27 27 27 ...
$ Lat_Min : num 44.3 53.5 43.6 58.3 58.3 ...
$ Lon : int 82 82 82 82 82 82 82 82 82 82 ...
$ Lon_Deg : int -82 -82 -82 -82 -82 -82 -82 -82 -82 -82 ...
$ Lon_Min : num 38.4 32 28.8 48.3 48.3 ...
$ Region.y : chr "TB" "MW" "MW" "MW" ...
$ Location.y : chr "big bayou" "found at gandy ramp" "Ruskin" "Clearwater" ...
$ LatDD : num 27.7 27.9 27.7 28 28 ...
$ LonDD : num -82.6 -82.5 -82.5 -82.8 -82.8 ...
$ INTSITE : num 770 770 770 770 770 770 770 770 770 770 ...
$ LonDD.1 : num 136935 147519 152596 120875 120875 ...
$ LatDD.1 : num 377392 394346 376032 403396 403396 ...
mercFL25miA <- DFfl25miA %>%
group_by(INTSITE,Species) %>%
summarise(n = n(), aveWt = mean(Whole_weight, na.rm = TRUE),
aveLt = mean(TL, na.rm = TRUE), aveHg = mean(Hg, na.rm = TRUE))
mercFL25miA
DFfl100miA$Whole_weight <- as.numeric(DFfl100miA$Whole_weight)
NAs introduced by coercion
DFfl100miA$TL <- as.numeric(DFfl100miA$TL)
mercFL100miA <- DFfl100miA %>%
group_by(INTSITE,Species) %>%
summarise(n = n(), aveWt = mean(Whole_weight, na.rm = TRUE),
aveLt = mean(TL, na.rm = TRUE), aveHg = mean(Hg, na.rm = TRUE))
mercFL100miA
DFfl25miB$Whole_weight <- as.numeric(DFfl25miB$Whole_weight)
NAs introduced by coercion
DFfl25miB$TL <- as.numeric(DFfl25miB$TL)
mercFL25miB <- DFfl25miB %>%
group_by(INTSITE,Species) %>%
summarise(n = n(), aveWt = mean(Whole_weight, na.rm = TRUE),
aveLt = mean(TL, na.rm = TRUE), aveHg = mean(Hg, na.rm = TRUE))
mercFL25miB
DFfl100miB$Whole_weight <- as.numeric(DFfl100miB$Whole_weight)
NAs introduced by coercion
DFfl100miB$TL <- as.numeric(DFfl100miB$TL)
mercFL100miB <- DFfl100miB %>%
group_by(INTSITE,Species) %>%
summarise(n = n(), aveWt = mean(Whole_weight, na.rm = TRUE),
aveLt = mean(TL, na.rm = TRUE), aveHg = mean(Hg, na.rm = TRUE))
mercFL100miB
#MRIP
landTtimA
landTtimB
mripFLa <- landTtimA %>%
group_by(INTSITE,common) %>%
summarise(n = n(), landings = sum(landing), biomass = sum(wgt_ab1, na.rm = TRUE))
mripFLa
mripFLb <- landTtimB %>%
group_by(INTSITE,common) %>%
summarise(n = n(), landings = sum(landing), biomass = sum(wgt_ab1, na.rm = TRUE))
mripFLb
Add in NOLA
str(DFla25miA)
'data.frame': 256 obs. of 24 variables:
$ WaterBody : chr "Grand Bayou Blue southwest of Galliano, Louisiana" "Grand Bayou Blue southwest of Galliano, Louisiana" "Grand Bayou Blue southwest of Galliano, Louisiana" "Grand Bayou Blue southwest of Galliano, Louisiana" ...
$ SiteNum : chr "3438" "3438" "3438" "3438" ...
$ Subsegment : chr "LA120706_00" "LA120706_00" "LA120706_00" "LA120706_00" ...
$ CollectDate : chr "2/5/2007" "2/5/2007" "2/5/2007" "2/5/2007" ...
$ CollectYear : int 2007 2007 2007 2007 2007 2007 2007 2007 2007 2007 ...
$ SampleType : chr "Single Fish" "Single Fish" "Composite" "Composite" ...
$ Parameter : chr "MERCURY" "MERCURY" "MERCURY" "MERCURY" ...
$ HgResult : num 0.0864 0.0692 0.0303 0.0316 0.0452 0.0233 0.095 0.0537 0.0843 0.0771 ...
$ Units : chr "ppm" "ppm" "ppm" "ppm" ...
$ spName : chr "SHEEPSHEAD" "SPOTTED SEATROUT" "SHEEPSHEAD" "SHEEPSHEAD" ...
$ aveWt_g : num 1590 385 594 918 1103 ...
$ aveLen_cm : num 45.4 34 32.1 37.4 38.8 38.5 50.7 40.1 54.7 62.1 ...
$ PROJ_NUM : chr "WQ1994001" "WQ1994001" "WQ1994001" "WQ1994001" ...
$ PROJECT_NAME : chr "Mercury Contaminant Levels in Louisiana Biota" "Mercury Contaminant Levels in Louisiana Biota" "Mercury Contaminant Levels in Louisiana Biota" "Mercury Contaminant Levels in Louisiana Biota" ...
$ SITE_NAME : chr "Grand Bayou Blue southwest of Galliano, Louisiana" "Grand Bayou Blue southwest of Galliano, Louisiana" "Grand Bayou Blue southwest of Galliano, Louisiana" "Grand Bayou Blue southwest of Galliano, Louisiana" ...
$ SITE_LOCATION: chr "4.6 miles southwest of Galliano, 5.3 miles northwest of Golden Meadow, 12.7 miles south of Larose" "4.6 miles southwest of Galliano, 5.3 miles northwest of Golden Meadow, 12.7 miles south of Larose" "4.6 miles southwest of Galliano, 5.3 miles northwest of Golden Meadow, 12.7 miles south of Larose" "4.6 miles southwest of Galliano, 5.3 miles northwest of Golden Meadow, 12.7 miles south of Larose" ...
$ SUBSEGMENT : chr "LA120706_00" "LA120706_00" "LA120706_00" "LA120706_00" ...
$ UTM_E : num 757506 757506 757506 757506 757506 ...
$ UTM_N : num 3254192 3254192 3254192 3254192 3254192 ...
$ LONG : num -90.3 -90.3 -90.3 -90.3 -90.3 ...
$ LAT : num 29.4 29.4 29.4 29.4 29.4 ...
$ INTSITE : num 222 222 222 222 222 222 222 222 222 222 ...
$ LONG.1 : num 1095774 1095774 1095774 1095774 1095774 ...
$ LAT.1 : num 99135 99135 99135 99135 99135 ...
mercLA25miA <- DFla25miA %>%
group_by(INTSITE,spName) %>%
summarise(n = n(), aveWt = mean(aveWt_g, na.rm = TRUE),
aveLt = mean(aveLen_cm, na.rm = TRUE), aveHg = mean(HgResult, na.rm = TRUE))
mercLA25miA
mercLA100miA <- DFla100miA %>%
group_by(INTSITE,spName) %>%
summarise(n = n(), aveWt = mean(aveWt_g, na.rm = TRUE),
aveLt = mean(aveLen_cm, na.rm = TRUE), aveHg = mean(HgResult, na.rm = TRUE))
mercLA100miA
mercLA25miBA <- DFla25miBA %>%
group_by(INTSITE,spName) %>%
summarise(n = n(), aveWt = mean(aveWt_g, na.rm = TRUE),
aveLt = mean(aveLen_cm, na.rm = TRUE), aveHg = mean(HgResult, na.rm = TRUE))
mercLA25miBA
mercLA100miBA <- DFla100miBA %>%
group_by(INTSITE,spName) %>%
summarise(n = n(), aveWt = mean(aveWt_g, na.rm = TRUE),
aveLt = mean(aveLen_cm, na.rm = TRUE), aveHg = mean(HgResult, na.rm = TRUE))
mercLA100miBA
#MRIP
mripLAa <- landOtimA %>%
group_by(INTSITE,common) %>%
summarise(n = n(), landings = sum(landing), biomass = sum(wgt_ab1, na.rm = TRUE))
mripLAa
mripLAb <- landOtimB %>%
group_by(INTSITE,common) %>%
summarise(n = n(), landings = sum(landing), biomass = sum(wgt_ab1, na.rm = TRUE))
mripLAb
NA
Try to indicate which entries are tested in mrip
str(landOtimA)
'data.frame': 545 obs. of 20 variables:
$ YEAR : chr "2005" "2005" "2005" "2005" ...
$ REG_RES : int 7 7 7 7 7 7 7 7 7 7 ...
$ ST_RES : int 22 22 22 22 22 22 22 22 22 22 ...
$ CNTY_RES : int 95 95 95 51 51 51 51 51 51 51 ...
$ ST : chr "22" "22" "22" "22" ...
$ CNTY : int 57 57 57 57 57 57 57 57 57 51 ...
$ INTSITE : int 231 231 231 231 231 231 231 231 231 222 ...
$ MODE_F : int 5 5 5 5 5 5 5 5 5 5 ...
$ MODE_FX : chr "3" "3" "3" "3" ...
$ ID_CODE : chr "1517620051011001" "1517620051011001" "1517620051011001" "1517620051011003" ...
$ ZIP : chr "70084" "70084" "70084" "70062" ...
$ common : chr "SHEEPSHEAD" "SPOTTED SEATROUT" "RED DRUM" "SPOTTED SEATROUT" ...
$ WAVE : int 5 5 5 5 5 5 5 5 5 5 ...
$ SP_CODE : num 8.84e+09 8.84e+09 8.84e+09 8.84e+09 8.84e+09 ...
$ harv_d3 : num 0 0 0 0 0 0 0 0 0 0 ...
$ wgt_ab1 : num 0.45 0.275 8.7 0.254 0.798 ...
$ landing : num 0.5 1 4.5 0.725 0.725 ...
$ testStatus : chr "Not Tested" "Not Tested" "Not Tested" "Not Tested" ...
$ testStatus25mi : chr "Tested" "Tested" "Tested" "Tested" ...
$ testStatus100mi: chr "Tested" "Tested" "Tested" "Tested" ...
str(mercLA25miA)
Classes ‘grouped_df’, ‘tbl_df’, ‘tbl’ and 'data.frame': 25 obs. of 6 variables:
$ INTSITE: num 151 151 151 151 151 155 155 155 155 155 ...
$ spName : chr "BLACK DRUM" "RED DRUM" "SHEEPSHEAD" "SOUTHERN FLOUNDER" ...
$ n : int 3 8 8 3 6 4 27 9 5 5 ...
$ aveWt : num 404 1756 993 467 445 ...
$ aveLt : num 193 53.2 37.1 33.6 35.8 ...
$ aveHg : num 0.0336 0.0753 0.054 0.0251 0.0667 ...
- attr(*, "groups")=Classes ‘tbl_df’, ‘tbl’ and 'data.frame': 5 obs. of 2 variables:
..$ INTSITE: num 151 155 222 231 306
..$ .rows :List of 5
.. ..$ : int 1 2 3 4 5
.. ..$ : int 6 7 8 9 10
.. ..$ : int 11 12 13 14 15
.. ..$ : int 16 17 18 19 20
.. ..$ : int 21 22 23 24 25
..- attr(*, ".drop")= logi TRUE
unique(mercLA25miA$spName)
[1] "BLACK DRUM" "RED DRUM" "SHEEPSHEAD" "SOUTHERN FLOUNDER"
[5] "SPOTTED SEATROUT"
#top sites = 222, 231, 306, 155, 151, within 25 mi
unique(mercLA25miA$INTSITE)
[1] 151 155 222 231 306
unique(landOtimA$INTSITE)
[1] 231 222 306 155 151
O222a25 <- subset(mercLA25miA,mercLA25miA$INTSITE == '222')
O231a25 <- subset(mercLA25miA,mercLA25miA$INTSITE == '231')
O306a25 <- subset(mercLA25miA,mercLA25miA$INTSITE == '306')
O155a25 <- subset(mercLA25miA,mercLA25miA$INTSITE == '155')
O151a25 <- subset(mercLA25miA,mercLA25miA$INTSITE == '151')
indLA25a222 <- which((landOtimA$common %in% O222a25$spName) & (landOtimA$INTSITE == '222'))
str(indLA25a222)
int [1:173] 11 12 19 22 23 24 31 37 38 40 ...
landOtimA[indLA25a222,]
landOtimA$testStatus25mi <- rep('NA', nrow(landOtimA))
landOtimA$testStatus25mi[indLA25a222] <- 'Tested'
indLA25a231 <- which((landOtimA$common %in% O231a25$spName) & (landOtimA$INTSITE == '231'))
str(indLA25a231)
int [1:101] 1 2 3 4 5 6 8 9 13 14 ...
landOtimA[indLA25a231,]
landOtimA$testStatus25mi[indLA25a231] <- 'Tested'
indLA25a306 <- which((landOtimA$common %in% O306a25$spName) & (landOtimA$INTSITE == '306'))
str(indLA25a306)
int [1:39] 58 175 176 177 178 191 192 194 215 216 ...
landOtimA[indLA25a306,]
landOtimA$testStatus25mi[indLA25a306] <- 'Tested'
indLA25a155 <- which((landOtimA$common %in% O155a25$spName) & (landOtimA$INTSITE == '155'))
str(indLA25a155)
int [1:30] 64 145 146 148 149 150 203 204 209 301 ...
landOtimA[indLA25a155,]
landOtimA$testStatus25mi[indLA25a155] <- 'Tested'
indLA25a151 <- which((landOtimA$common %in% O151a25$spName) & (landOtimA$INTSITE == '151'))
str(indLA25a151)
int [1:11] 236 250 291 292 322 410 411 519 525 537 ...
landOtimA[indLA25a151,]
landOtimA$testStatus25mi[indLA25a151] <- 'Tested'
landOtimA$testStatus25mi[landOtimA$testStatus25mi %in% 'NA'] <- 'Not Tested'
#top sites = 222, 231, 306, 155, 151, within 100 mi
unique(mercLA100miA$INTSITE)
[1] 151 155 222 231 306
unique(mercLA100miA$spName)
[1] "BLACK DRUM" "BLACKFIN TUNA" "COBIA" "GRAY SNAPPER"
[5] "GRAY TRIGGERFISH" "KING MACKEREL" "RED DRUM" "RED SNAPPER"
[9] "SHEEPSHEAD" "SOUTHERN FLOUNDER" "SPOTTED SEATROUT" "STRIPED BASS"
unique(landOtimA$INTSITE)
[1] 231 222 306 155 151
O222a100 <- subset(mercLA100miA,mercLA100miA$INTSITE == '222')
O231a100 <- subset(mercLA100miA,mercLA100miA$INTSITE == '231')
O306a100 <- subset(mercLA100miA,mercLA100miA$INTSITE == '306')
O155a100 <- subset(mercLA100miA,mercLA100miA$INTSITE == '155')
O151a100 <- subset(mercLA100miA,mercLA100miA$INTSITE == '151')
indLA100a222 <- which((landOtimA$common %in% O222a100$spName) & (landOtimA$INTSITE == '222'))
str(indLA100a222)
int [1:175] 11 12 19 22 23 24 31 37 38 40 ...
landOtimA[indLA100a222,]
landOtimA$testStatus100mi <- rep('NA', nrow(landOtimA))
landOtimA$testStatus100mi[indLA100a222] <- 'Tested'
indLA100a231 <- which((landOtimA$common %in% O231a100$spName) & (landOtimA$INTSITE == '231'))
str(indLA100a231)
int [1:101] 1 2 3 4 5 6 8 9 13 14 ...
landOtimA[indLA100a231,]
landOtimA$testStatus100mi[indLA100a231] <- 'Tested'
indLA100a306 <- which((landOtimA$common %in% O306a100$spName) & (landOtimA$INTSITE == '306'))
str(indLA100a306)
int [1:39] 58 175 176 177 178 191 192 194 215 216 ...
landOtimA[indLA100a306,]
landOtimA$testStatus100mi[indLA100a306] <- 'Tested'
indLA100a155 <- which((landOtimA$common %in% O155a100$spName) & (landOtimA$INTSITE == '155'))
str(indLA100a155)
int [1:30] 64 145 146 148 149 150 203 204 209 301 ...
landOtimA[indLA100a155,]
landOtimA$testStatus100mi[indLA100a155] <- 'Tested'
indLA100a151 <- which((landOtimA$common %in% O151a100$spName) & (landOtimA$INTSITE == '151'))
str(indLA100a151)
int [1:11] 236 250 291 292 322 410 411 519 525 537 ...
landOtimA[indLA100a151,]
landOtimA$testStatus100mi[indLA100a151] <- 'Tested'
landOtimA$testStatus100mi[landOtimA$testStatus100mi %in% 'NA'] <- 'Not Tested'
landOtimA
LAaProp25mi <- landOtimA %>%
group_by(INTSITE, testStatus25mi) %>%
summarise(n_distinct(common), landings = sum(landing), biomass = sum(wgt_ab1))
LAaProp100mi <- landOtimA %>%
group_by(INTSITE, testStatus100mi) %>%
summarise(n_distinct(common), landings = sum(landing), biomass = sum(wgt_ab1))
write.csv(LAaProp25mi, "~/FUIteam/PydioData/env/data_outputs/LAaProp25mi.csv")
write.csv(LAaProp100mi, "~/FUIteam/PydioData/env/data_outputs/LAaProp100mi.csv")
New Orleans Time Chunk B
str(landOtimB)
'data.frame': 225 obs. of 19 variables:
$ YEAR : chr "2010" "2010" "2010" "2010" ...
$ REG_RES : int 7 7 7 7 7 7 7 7 7 7 ...
$ ST_RES : int 22 22 22 22 22 22 22 22 22 22 ...
$ CNTY_RES : int 89 51 93 51 51 51 89 71 71 71 ...
$ ST : chr "22" "22" "22" "22" ...
$ CNTY : int 51 51 51 51 51 51 51 51 51 51 ...
$ INTSITE : int 159 222 222 222 222 222 222 222 222 222 ...
$ MODE_F : int 5 5 5 5 5 5 3 5 5 5 ...
$ MODE_FX : chr "3" "3" "3" "3" ...
$ ID_CODE : chr "1541220100214013" "1517620100130002" "1517620100130005" "1517620100214009" ...
$ ZIP : chr "70070" "70358" "70090" "70358" ...
$ common : chr "RED DRUM" "SPOTTED SEATROUT" "SPOTTED SEATROUT" "SHEEPSHEAD" ...
$ WAVE : int 1 1 1 1 1 1 2 2 2 2 ...
$ SP_CODE : num 8.84e+09 8.84e+09 8.84e+09 8.84e+09 8.84e+09 ...
$ harv_d3 : num 0 0 9 2.65 0 ...
$ wgt_ab1 : num 1.8 5.56 2.88 2.12 3.4 ...
$ landing : num 1 15.89 9 2.65 9.5 ...
$ testStatus25mi : chr "Tested" "Tested" "Tested" "Tested" ...
$ testStatus100mi: chr "Tested" "Tested" "Tested" "Tested" ...
str(mercLA25miBA)
Classes ‘grouped_df’, ‘tbl_df’, ‘tbl’ and 'data.frame': 23 obs. of 6 variables:
$ INTSITE: chr "151ba" "151ba" "151ba" "151ba" ...
$ spName : chr "BLACK DRUM" "RED DRUM" "SHEEPSHEAD" "SOUTHERN FLOUNDER" ...
$ n : int 3 8 8 3 6 5 12 11 6 10 ...
$ aveWt : num 404 1756 993 467 445 ...
$ aveLt : num 193 53.2 37.1 33.6 35.8 ...
$ aveHg : num 0.0336 0.0753 0.054 0.0251 0.0667 ...
- attr(*, "groups")=Classes ‘tbl_df’, ‘tbl’ and 'data.frame': 5 obs. of 2 variables:
..$ INTSITE: chr "151ba" "159ba" "222ba" "306ba" ...
..$ .rows :List of 5
.. ..$ : int 1 2 3 4 5
.. ..$ : int 6 7 8 9 10
.. ..$ : int 11 12 13 14 15
.. ..$ : int 16 17 18 19 20
.. ..$ : int 21 22 23
..- attr(*, ".drop")= logi TRUE
unique(mercLA25miBA$spName)
[1] "BLACK DRUM" "RED DRUM" "SHEEPSHEAD" "SOUTHERN FLOUNDER"
[5] "SPOTTED SEATROUT" "STRIPED BASS"
str(mercLA25miA)
Classes ‘grouped_df’, ‘tbl_df’, ‘tbl’ and 'data.frame': 25 obs. of 6 variables:
$ INTSITE: num 151 151 151 151 151 155 155 155 155 155 ...
$ spName : chr "BLACK DRUM" "RED DRUM" "SHEEPSHEAD" "SOUTHERN FLOUNDER" ...
$ n : int 3 8 8 3 6 4 27 9 5 5 ...
$ aveWt : num 404 1756 993 467 445 ...
$ aveLt : num 193 53.2 37.1 33.6 35.8 ...
$ aveHg : num 0.0336 0.0753 0.054 0.0251 0.0667 ...
- attr(*, "groups")=Classes ‘tbl_df’, ‘tbl’ and 'data.frame': 5 obs. of 2 variables:
..$ INTSITE: num 151 155 222 231 306
..$ .rows :List of 5
.. ..$ : int 1 2 3 4 5
.. ..$ : int 6 7 8 9 10
.. ..$ : int 11 12 13 14 15
.. ..$ : int 16 17 18 19 20
.. ..$ : int 21 22 23 24 25
..- attr(*, ".drop")= logi TRUE
str(mercLA100miA)
Classes ‘grouped_df’, ‘tbl_df’, ‘tbl’ and 'data.frame': 58 obs. of 6 variables:
$ INTSITE: num 151 151 151 151 151 151 151 151 151 151 ...
$ spName : chr "BLACK DRUM" "BLACKFIN TUNA" "COBIA" "GRAY SNAPPER" ...
$ n : int 14 1 3 17 2 5 145 3 73 28 ...
$ aveWt : num 950 9049 13922 1930 1108 ...
$ aveLt : num 73.2 77.5 113.4 49 40.6 ...
$ aveHg : num 0.0794 0.7111 0.5897 0.1207 0.1118 ...
- attr(*, "groups")=Classes ‘tbl_df’, ‘tbl’ and 'data.frame': 5 obs. of 2 variables:
..$ INTSITE: num 151 155 222 231 306
..$ .rows :List of 5
.. ..$ : int 1 2 3 4 5 6 7 8 9 10 ...
.. ..$ : int 13 14 15 16 17 18 19 20 21 22 ...
.. ..$ : int 24 25 26 27 28 29 30 31 32 33 ...
.. ..$ : int 36 37 38 39 40 41 42 43 44 45 ...
.. ..$ : int 48 49 50 51 52 53 54 55 56 57 ...
..- attr(*, ".drop")= logi TRUE
unique(mercLA25miBA$spName)
[1] "BLACK DRUM" "RED DRUM" "SHEEPSHEAD" "SOUTHERN FLOUNDER"
[5] "SPOTTED SEATROUT" "STRIPED BASS"
unique(mercLA100miBA$spName)
[1] "BLACK DRUM" "BLACKFIN TUNA" "COBIA" "GRAY SNAPPER"
[5] "GRAY TRIGGERFISH" "KING MACKEREL" "RED DRUM" "RED SNAPPER"
[9] "SHEEPSHEAD" "SOUTHERN FLOUNDER" "SPOTTED SEATROUT" "STRIPED BASS"
unique(landOtimB$common)
[1] "RED DRUM" "SPOTTED SEATROUT" "SHEEPSHEAD"
[4] "SOUTHERN FLOUNDER" "HARDHEAD CATFISH" "GAFFTOPSAIL CATFISH"
[7] "SAND SEATROUT" "PINFISH" "SOUTHERN KINGFISH"
[10] "BLACK DRUM" "ATLANTIC CROAKER" "GULF KINGFISH"
[13] "SILVER PERCH" "BLUE CATFISH" "BLUEGILL"
[16] "SKIPJACK HERRING"
#top sites = 222, 159, 306, 3325, 151, within 25 mi
unique(mercLA25miBA$INTSITE)
[1] "151ba" "159ba" "222ba" "306ba" "3325ba"
unique(landOtimB$INTSITE)
[1] 159 222 3325 151 306
O222ba25 <- subset(mercLA25miBA,mercLA25miBA$INTSITE == '222ba')
O3325ba25 <- subset(mercLA25miBA,mercLA25miBA$INTSITE == '3325ba')
O306ba25 <- subset(mercLA25miBA,mercLA25miBA$INTSITE == '306ba')
O159ba25 <- subset(mercLA25miBA,mercLA25miBA$INTSITE == '159ba')
O151ba25 <- subset(mercLA25miBA,mercLA25miBA$INTSITE == '151ba')
indLA25b222 <- which((landOtimB$common %in% O222ba25$spName) & (landOtimB$INTSITE == '222'))
str(indLA25b222)
int [1:116] 2 3 4 5 6 7 8 9 10 21 ...
landOtimB[indLA25b222,]
landOtimB$testStatus25mi <- rep('NA', nrow(landOtimB))
landOtimB$testStatus25mi[indLA25b222] <- 'Tested'
subset(landOtimB, landOtimB$INTSITE == '3325')
indLA25b3325 <- which((landOtimB$common %in% O3325ba25$spName) & (landOtimB$INTSITE == '3325'))
str(indLA25b3325)
int [1:4] 12 13 122 217
landOtimB[indLA25b3325,]
landOtimB$testStatus25mi[indLA25b3325] <- 'Tested'
indLA25b306 <- which((landOtimB$common %in% O306ba25$spName) & (landOtimB$INTSITE == '306'))
str(indLA25b306)
int [1:4] 88 89 152 153
landOtimB[indLA25b306,]
landOtimB$testStatus25mi[indLA25b306] <- 'Tested'
indLA25b159 <- which((landOtimB$common %in% O159ba25$spName) & (landOtimB$INTSITE == '159'))
str(indLA25b159)
int [1:19] 1 66 75 77 79 81 90 99 100 101 ...
landOtimB[indLA25b159,]
landOtimB$testStatus25mi[indLA25b159] <- 'Tested'
indLA25b151 <- which((landOtimB$common %in% O151ba25$spName) & (landOtimB$INTSITE == '151'))
str(indLA25b151)
int [1:6] 26 120 159 160 177 181
landOtimB[indLA25b151,]
landOtimB$testStatus25mi[indLA25b151] <- 'Tested'
landOtimB$testStatus25mi[landOtimB$testStatus25mi %in% 'NA'] <- 'Not Tested'
#top sites = 222, 159, 306, 3325, 151, within 100 mi
unique(mercLA100miBA$INTSITE)
[1] "151ba" "159ba" "222ba" "306ba" "3325ba"
str(mercLA25miBA)
Classes ‘grouped_df’, ‘tbl_df’, ‘tbl’ and 'data.frame': 23 obs. of 6 variables:
$ INTSITE: chr "151ba" "151ba" "151ba" "151ba" ...
$ spName : chr "BLACK DRUM" "RED DRUM" "SHEEPSHEAD" "SOUTHERN FLOUNDER" ...
$ n : int 3 8 8 3 6 5 12 11 6 10 ...
$ aveWt : num 404 1756 993 467 445 ...
$ aveLt : num 193 53.2 37.1 33.6 35.8 ...
$ aveHg : num 0.0336 0.0753 0.054 0.0251 0.0667 ...
- attr(*, "groups")=Classes ‘tbl_df’, ‘tbl’ and 'data.frame': 5 obs. of 2 variables:
..$ INTSITE: chr "151ba" "159ba" "222ba" "306ba" ...
..$ .rows :List of 5
.. ..$ : int 1 2 3 4 5
.. ..$ : int 6 7 8 9 10
.. ..$ : int 11 12 13 14 15
.. ..$ : int 16 17 18 19 20
.. ..$ : int 21 22 23
..- attr(*, ".drop")= logi TRUE
unique(mercLA100miBA$spName)
[1] "BLACK DRUM" "BLACKFIN TUNA" "COBIA" "GRAY SNAPPER"
[5] "GRAY TRIGGERFISH" "KING MACKEREL" "RED DRUM" "RED SNAPPER"
[9] "SHEEPSHEAD" "SOUTHERN FLOUNDER" "SPOTTED SEATROUT" "STRIPED BASS"
unique(landOtimB$INTSITE)
[1] 159 222 3325 151 306
O222ba100 <- subset(mercLA100miBA,mercLA100miBA$INTSITE == '222ba')
O3325ba100 <- subset(mercLA100miBA,mercLA100miBA$INTSITE == '3325ba')
O306ba100 <- subset(mercLA100miBA,mercLA100miBA$INTSITE == '306ba')
O159ba100 <- subset(mercLA100miBA,mercLA100miBA$INTSITE == '159ba')
O151ba100 <- subset(mercLA100miBA,mercLA100miBA$INTSITE == '151ba')
indLA100b222 <- which((landOtimB$common %in% O222ba100$spName) & (landOtimB$INTSITE == '222'))
str(indLA100b222)
int [1:116] 2 3 4 5 6 7 8 9 10 21 ...
landOtimB[indLA100b222,]
landOtimB$testStatus100mi <- rep('NA', nrow(landOtimB))
landOtimB$testStatus100mi[indLA100b222] <- 'Tested'
subset(landOtimB, landOtimB$INTSITE == '3325')
indLA100b3325 <- which((landOtimB$common %in% O3325ba100$spName) & (landOtimB$INTSITE == '3325'))
str(indLA100b3325)
int [1:4] 12 13 122 217
landOtimB[indLA100b3325,]
landOtimB$testStatus100mi[indLA100b3325] <- 'Tested'
indLA100b306 <- which((landOtimB$common %in% O306ba100$spName) & (landOtimB$INTSITE == '306'))
str(indLA100b306)
int [1:4] 88 89 152 153
landOtimB[indLA100b306,]
landOtimB$testStatus100mi[indLA100b306] <- 'Tested'
indLA100b159 <- which((landOtimB$common %in% O159ba100$spName) & (landOtimB$INTSITE == '159'))
str(indLA100b159)
int [1:19] 1 66 75 77 79 81 90 99 100 101 ...
landOtimB[indLA100b159,]
landOtimB$testStatus100mi[indLA100b159] <- 'Tested'
indLA100b151 <- which((landOtimB$common %in% O151ba100$spName) & (landOtimB$INTSITE == '151'))
str(indLA100b151)
int [1:6] 26 120 159 160 177 181
landOtimB[indLA100b151,]
landOtimB$testStatus100mi[indLA100b151] <- 'Tested'
landOtimB$testStatus100mi[landOtimB$testStatus100mi %in% 'NA'] <- 'Not Tested'
unique(landOtimB$testStatus100mi)
[1] "Tested" "Not Tested"
LAbaProp25mi <- landOtimB %>%
group_by(INTSITE, testStatus25mi) %>%
summarise(n_distinct(common), landings = sum(landing), biomass = sum(wgt_ab1))
LAbaProp100mi <- landOtimB %>%
group_by(INTSITE, testStatus100mi) %>%
summarise(n_distinct(common), landings = sum(landing), biomass = sum(wgt_ab1))
write.csv(LAbaProp25mi, "~/FUIteam/PydioData/env/data_outputs/LAbaProp25mi.csv")
write.csv(LAbaProp100mi, "~/FUIteam/PydioData/env/data_outputs/LAbaProp100mi.csv")
Now onto Florida
str(landTtimA)
'data.frame': 1622 obs. of 19 variables:
$ YEAR : chr "2005" "2005" "2005" "2005" ...
$ REG_RES : int 7 7 7 7 7 7 7 7 7 7 ...
$ ST_RES : int 12 12 12 12 12 12 12 12 12 12 ...
$ CNTY_RES : int 103 103 103 103 103 103 103 103 103 57 ...
$ ST : chr "12" "12" "12" "12" ...
$ CNTY : int 103 103 103 103 103 103 103 103 103 103 ...
$ INTSITE : int 615 770 615 615 615 770 770 769 769 769 ...
$ MODE_F : int 1 1 1 1 1 1 1 1 1 1 ...
$ MODE_FX : chr "3" "3" "3" "3" ...
$ ID_CODE : chr "1297920050326028" "1520920050324001" "1520920050317005" "1520920050317005" ...
$ ZIP : chr "33710" "33706" "33708" "33708" ...
$ common : chr "SPANISH MACKEREL" "SCALED SARDINE" "SPOTTED SEATROUT" "SILVER PERCH" ...
$ WAVE : int 2 2 2 2 2 1 1 1 1 1 ...
$ SP_CODE : num 8.85e+09 8.75e+09 8.84e+09 8.84e+09 8.84e+09 ...
$ harv_d3 : num 0 0 3 10 0 0 0 0 0 0 ...
$ wgt_ab1 : num 3.779 0.413 2.86 1.551 0.602 ...
$ landing : num 9.04 24 3 18.03 2.01 ...
$ testStatus25mi : chr "Tested" "Not Tested" "Not Tested" "Not Tested" ...
$ testStatus100mi: chr "Tested" "Not Tested" "Not Tested" "Not Tested" ...
str(mercFL25miA)
Classes ‘grouped_df’, ‘tbl_df’, ‘tbl’ and 'data.frame': 67 obs. of 6 variables:
$ INTSITE: num 615 615 615 615 615 615 615 615 615 615 ...
$ Species: chr "BULL SHARK" "COBIA" "GAG" "GRAY SNAPPER" ...
$ n : int 1 1 2 9 8 8 1 5 1 1 ...
$ aveWt : num NaN NaN 6 2 1 5.75 NaN 2.2 2 3 ...
$ aveLt : num 2800 1279 714 528 393 ...
$ aveHg : num 1.9 0.8 0.48 0.323 0.124 ...
- attr(*, "groups")=Classes ‘tbl_df’, ‘tbl’ and 'data.frame': 5 obs. of 2 variables:
..$ INTSITE: num 615 632 742 769 770
..$ .rows :List of 5
.. ..$ : int 1 2 3 4 5 6 7 8 9 10 ...
.. ..$ : int 14 15 16 17 18 19 20 21 22 23 ...
.. ..$ : int 26 27 28 29 30 31 32 33 34 35 ...
.. ..$ : int 38 39 40 41 42 43 44 45 46 47 ...
.. ..$ : int 53 54 55 56 57 58 59 60 61 62 ...
..- attr(*, ".drop")= logi TRUE
str(mercFL100miA)
Classes ‘grouped_df’, ‘tbl_df’, ‘tbl’ and 'data.frame': 110 obs. of 6 variables:
$ INTSITE: num 615 615 615 615 615 615 615 615 615 615 ...
$ Species: chr "" "ATLANTIC SHARPNOSE SHARK" "BLACKTIP SHARK" "BONNETHEAD" ...
$ n : int 1 1 3 11 5 3 11 11 3 9 ...
$ aveWt : num NaN 117 2530 3581 4648 ...
$ aveLt : num 3429 310 1115 916 2103 ...
$ aveHg : num 3.5 0.088 1.463 0.789 1.978 ...
- attr(*, "groups")=Classes ‘tbl_df’, ‘tbl’ and 'data.frame': 5 obs. of 2 variables:
..$ INTSITE: num 615 632 742 769 770
..$ .rows :List of 5
.. ..$ : int 1 2 3 4 5 6 7 8 9 10 ...
.. ..$ : int 24 25 26 27 28 29 30 31 32 33 ...
.. ..$ : int 45 46 47 48 49 50 51 52 53 54 ...
.. ..$ : int 69 70 71 72 73 74 75 76 77 78 ...
.. ..$ : int 90 91 92 93 94 95 96 97 98 99 ...
..- attr(*, ".drop")= logi TRUE
unique(mercFL25miA$Species)
[1] "BULL SHARK" "COBIA" "GAG" "GRAY SNAPPER"
[5] "GRAY TRIGGERFISH" "GREATER AMBERJACK" "GULF FLOUNDER" "KING MACKEREL"
[9] "LITTLE TUNNY" "RED DRUM" "RED GROUPER" "SNOOK"
[13] "SPANISH MACKEREL" "LANE SNAPPER" "VERMILION SNAPPER"
unique(mercFL100miA$Species)
[1] "" "ATLANTIC SHARPNOSE SHARK" "BLACKTIP SHARK"
[4] "BONNETHEAD" "BULL SHARK" "COBIA"
[7] "CREVALLE JACK" "FLORIDA POMPANO" "GAG"
[10] "GRAY SNAPPER" "GRAY TRIGGERFISH" "GREATER AMBERJACK"
[13] "GULF FLOUNDER" "KING MACKEREL" "LANE SNAPPER"
[16] "LEMON SHARK" "LITTLE TUNNY" "RED DRUM"
[19] "RED GROUPER" "SNOOK" "SOUTHERN KINGFISH"
[22] "SPANISH MACKEREL" "VERMILION SNAPPER" "SPOTTED SEATROUT"
#need to convert lower case to upper case for the two dataframes to be able search properly
mercFL25miA$Species <- toupper(mercFL25miA$Species)
mercFL25miA$Species <- trimws(mercFL25miA$Species, which = c("right"), whitespace = "[ \t\r\n]")
mercFL100miA$Species <- toupper(mercFL100miA$Species)
mercFL100miA$Species <- trimws(mercFL100miA$Species, which = c("right"), whitespace = "[ \t\r\n]")
#edit COMMON SNOOK in MRIP data to read SNOOK so they match
landTtimA$common <- gsub("COMMON SNOOK", "SNOOK", landTtimA$common)
#top sites = 770, 769, 632, 615, 742, within 25 mi
unique(mercFL25miA$INTSITE)
[1] 615 632 742 769 770
unique(landTtimA$INTSITE)
[1] 615 770 769 632 742
T615a25 <- subset(mercFL25miA,mercFL25miA$INTSITE == '615')
T632a25 <- subset(mercFL25miA,mercFL25miA$INTSITE == '632')
T742a25 <- subset(mercFL25miA,mercFL25miA$INTSITE == '742')
T769a25 <- subset(mercFL25miA,mercFL25miA$INTSITE == '769')
T770a25 <- subset(mercFL25miA,mercFL25miA$INTSITE == '770')
indFL25a615 <- which((landTtimA$common %in% T615a25$Species) & (landTtimA$INTSITE == '615'))
str(indFL25a615)
int [1:42] 1 121 124 149 274 276 277 398 421 885 ...
landTtimA[indFL25a615,]
unique((subset(landTtimA,landTtimA$INTSITE == '615'))$common)
[1] "SPANISH MACKEREL" "SPOTTED SEATROUT" "SILVER PERCH"
[4] "ATLANTIC CROAKER" "FLORIDA POMPANO" "GULF KINGFISH"
[7] "GRAY SNAPPER" "SILVER SEATROUT" "SOUTHERN KINGFISH"
[10] "SCALED SARDINE" "ATLANTIC BUMPER" "LADYFISH"
[13] "ATLANTIC THREAD HERRING" "SPANISH SARDINE" "PINFISH"
[16] "BLUEFISH" "MACKEREL GENUS" "HOGFISH"
[19] "ATLANTIC SHARPNOSE SHARK" "BLUE RUNNER" "WHITE GRUNT"
[22] "GULF FLOUNDER" "INSHORE LIZARDFISH" "BLACK DRUM"
[25] "SAND SEATROUT" "PIGFISH" "SPOTTAIL PINFISH"
[28] "SNOOK" "CREVALLE JACK" "LEFTEYE FLOUNDER FAMILY"
[31] "KINGFISH GENUS" "RED DRUM" "SEATROUT GENUS"
[34] "KING MACKEREL"
unique(T615a25$Species)
[1] "BULL SHARK" "COBIA" "GAG" "GRAY SNAPPER"
[5] "GRAY TRIGGERFISH" "GREATER AMBERJACK" "GULF FLOUNDER" "KING MACKEREL"
[9] "LITTLE TUNNY" "RED DRUM" "RED GROUPER" "SNOOK"
[13] "SPANISH MACKEREL"
landTtimA$testStatus25mi <- rep('NA', nrow(landTtimA))
landTtimA$testStatus25mi[indFL25a615] <- 'Tested'
indFL25a632 <- which((landTtimA$common %in% T632a25$Species) & (landTtimA$INTSITE == '632'))
str(indFL25a632)
int [1:59] 75 134 137 139 192 348 349 358 359 360 ...
landTtimA[indFL25a632,]
unique(T632a25$Species)
[1] "BULL SHARK" "COBIA" "GAG" "GRAY SNAPPER"
[5] "GRAY TRIGGERFISH" "GREATER AMBERJACK" "GULF FLOUNDER" "LANE SNAPPER"
[9] "RED DRUM" "SNOOK" "SPANISH MACKEREL" "VERMILION SNAPPER"
landTtimA$testStatus25mi[indFL25a632] <- 'Tested'
indFL25a742 <- which((landTtimA$common %in% T742a25$Species) & (landTtimA$INTSITE == '742'))
str(indFL25a742)
int [1:15] 54 189 196 198 214 308 369 370 516 1048 ...
landTtimA[indFL25a742,]
landTtimA$testStatus25mi[indFL25a742] <- 'Tested'
indFL25a769 <- which((landTtimA$common %in% T769a25$Species) & (landTtimA$INTSITE == '769'))
str(indFL25a769)
int [1:110] 111 112 113 114 115 116 117 118 119 120 ...
landTtimA[indFL25a769,]
landTtimA$testStatus25mi[indFL25a769] <- 'Tested'
indFL25a770 <- which((landTtimA$common %in% T770a25$Species) & (landTtimA$INTSITE == '770'))
str(indFL25a770)
int [1:129] 21 25 30 32 33 34 39 40 41 42 ...
landTtimA[indFL25a770,]
landTtimA$testStatus25mi[indFL25a770] <- 'Tested'
landTtimA$testStatus25mi[landTtimA$testStatus25mi %in% 'NA'] <- 'Not Tested'
100 mi chunk for FL, time frame A
#top sites = 770, 769, 632, 615, 742, within 100 mi
unique(mercFL100miA$INTSITE)
[1] 615 632 742 769 770
unique(landTtimA$INTSITE)
[1] 615 770 769 632 742
T615a100 <- subset(mercFL100miA,mercFL100miA$INTSITE == '615')
T632a100 <- subset(mercFL100miA,mercFL100miA$INTSITE == '632')
T742a100 <- subset(mercFL100miA,mercFL100miA$INTSITE == '742')
T769a100 <- subset(mercFL100miA,mercFL100miA$INTSITE == '769')
T770a100 <- subset(mercFL100miA,mercFL100miA$INTSITE == '770')
indFL100a615 <- which((landTtimA$common %in% T615a100$Species) & (landTtimA$INTSITE == '615'))
str(indFL100a615)
int [1:56] 1 121 122 124 149 258 261 274 276 277 ...
landTtimA[indFL100a615,]
unique((subset(landTtimA,landTtimA$INTSITE == '615'))$common)
[1] "SPANISH MACKEREL" "SPOTTED SEATROUT" "SILVER PERCH"
[4] "ATLANTIC CROAKER" "FLORIDA POMPANO" "GULF KINGFISH"
[7] "GRAY SNAPPER" "SILVER SEATROUT" "SOUTHERN KINGFISH"
[10] "SCALED SARDINE" "ATLANTIC BUMPER" "LADYFISH"
[13] "ATLANTIC THREAD HERRING" "SPANISH SARDINE" "PINFISH"
[16] "BLUEFISH" "MACKEREL GENUS" "HOGFISH"
[19] "ATLANTIC SHARPNOSE SHARK" "BLUE RUNNER" "WHITE GRUNT"
[22] "GULF FLOUNDER" "INSHORE LIZARDFISH" "BLACK DRUM"
[25] "SAND SEATROUT" "PIGFISH" "SPOTTAIL PINFISH"
[28] "SNOOK" "CREVALLE JACK" "LEFTEYE FLOUNDER FAMILY"
[31] "KINGFISH GENUS" "RED DRUM" "SEATROUT GENUS"
[34] "KING MACKEREL"
unique(T615a100$Species)
[1] "" "ATLANTIC SHARPNOSE SHARK" "BLACKTIP SHARK"
[4] "BONNETHEAD" "BULL SHARK" "COBIA"
[7] "CREVALLE JACK" "FLORIDA POMPANO" "GAG"
[10] "GRAY SNAPPER" "GRAY TRIGGERFISH" "GREATER AMBERJACK"
[13] "GULF FLOUNDER" "KING MACKEREL" "LANE SNAPPER"
[16] "LEMON SHARK" "LITTLE TUNNY" "RED DRUM"
[19] "RED GROUPER" "SNOOK" "SOUTHERN KINGFISH"
[22] "SPANISH MACKEREL" "VERMILION SNAPPER"
landTtimA$testStatus100mi <- rep('NA', nrow(landTtimA))
landTtimA$testStatus100mi[indFL100a615] <- 'Tested'
indFL100a632 <- which((landTtimA$common %in% T632a100$Species) & (landTtimA$INTSITE == '632'))
str(indFL100a632)
int [1:72] 75 133 134 136 137 139 192 346 347 348 ...
landTtimA[indFL100a632,]
unique(T632a100$Species)
[1] "" "BLACKTIP SHARK" "BONNETHEAD" "BULL SHARK"
[5] "COBIA" "FLORIDA POMPANO" "GAG" "GRAY SNAPPER"
[9] "GRAY TRIGGERFISH" "GREATER AMBERJACK" "GULF FLOUNDER" "KING MACKEREL"
[13] "LANE SNAPPER" "LEMON SHARK" "LITTLE TUNNY" "RED DRUM"
[17] "RED GROUPER" "SNOOK" "SOUTHERN KINGFISH" "SPANISH MACKEREL"
[21] "VERMILION SNAPPER"
landTtimA$testStatus100mi[indFL100a632] <- 'Tested'
indFL100a742 <- which((landTtimA$common %in% T742a100$Species) & (landTtimA$INTSITE == '742'))
str(indFL100a742)
int [1:25] 24 54 126 128 185 189 194 196 198 214 ...
landTtimA[indFL100a742,]
landTtimA$testStatus100mi[indFL100a742] <- 'Tested'
indFL100a769 <- which((landTtimA$common %in% T769a100$Species) & (landTtimA$INTSITE == '769'))
str(indFL100a769)
int [1:116] 12 111 112 113 114 115 116 117 118 119 ...
landTtimA[indFL100a769,]
landTtimA$testStatus100mi[indFL100a769] <- 'Tested'
indFL100a770 <- which((landTtimA$common %in% T770a100$Species) & (landTtimA$INTSITE == '770'))
str(indFL100a770)
int [1:134] 21 25 30 32 33 34 39 40 41 42 ...
landTtimA[indFL100a770,]
landTtimA$testStatus100mi[indFL100a770] <- 'Tested'
landTtimA$testStatus100mi[landTtimA$testStatus100mi %in% 'NA'] <- 'Not Tested'
FLaProp25mi <- landTtimA %>%
group_by(INTSITE, testStatus25mi) %>%
summarise(n_distinct(common), landings = sum(landing), biomass = sum(wgt_ab1))
FLaProp100mi <- landTtimA %>%
group_by(INTSITE, testStatus100mi) %>%
summarise(n_distinct(common), landings = sum(landing), biomass = sum(wgt_ab1))
write.csv(FLaProp25mi, "~/FUIteam/PydioData/env/data_outputs/FLaProp25mi.csv")
write.csv(FLaProp100mi, "~/FUIteam/PydioData/env/data_outputs/FLaProp100mi.csv")
Timeframe B, 25 mi
str(landTtimB)
'data.frame': 956 obs. of 19 variables:
$ YEAR : chr "2011" "2011" "2011" "2011" ...
$ REG_RES : int 7 7 7 7 7 7 7 7 7 7 ...
$ ST_RES : int 12 12 12 12 12 12 12 12 12 12 ...
$ CNTY_RES : int 103 57 57 57 57 57 57 57 57 57 ...
$ ST : chr "12" "12" "12" "12" ...
$ CNTY : int 103 103 103 103 103 103 103 103 103 103 ...
$ INTSITE : int 770 769 769 769 769 769 769 769 769 769 ...
$ MODE_F : int 1 5 5 1 1 1 1 1 1 1 ...
$ MODE_FX : chr "3" "3" "3" "3" ...
$ ID_CODE : chr "1544320110218005" "1000920110227003" "1000920110227003" "1000920110227010" ...
$ ZIP : chr "33778" "33547" "33547" "33612" ...
$ common : chr "SPANISH MACKEREL" "SAND SEATROUT" "SPANISH MACKEREL" "SCALED SARDINE" ...
$ WAVE : int 1 1 1 1 1 1 1 1 1 1 ...
$ SP_CODE : num 8.85e+09 8.84e+09 8.85e+09 8.75e+09 8.84e+09 ...
$ harv_d3 : num 1 0 0 0 0 0 0 0 0 0 ...
$ wgt_ab1 : num 0.38 0.239 0.545 0.149 0.16 ...
$ landing : num 1 0.556 1.112 12 2 ...
$ testStatus25mi : chr "Tested" "Not Tested" "Tested" "Not Tested" ...
$ testStatus100mi: chr "Tested" "Not Tested" "Tested" "Not Tested" ...
str(mercFL25miB)
Classes ‘grouped_df’, ‘tbl_df’, ‘tbl’ and 'data.frame': 70 obs. of 6 variables:
$ INTSITE: chr "3802b" "3802b" "3802b" "3802b" ...
$ Species: chr "Bull shark " "Cobia " "Gag " "Gray snapper " ...
$ n : int 1 2 3 9 8 12 1 5 1 1 ...
$ aveWt : num NaN 7 4.5 2 1 ...
$ aveLt : num 2800 1120 687 528 393 ...
$ aveHg : num 1.9 0.72 0.493 0.323 0.124 ...
- attr(*, "groups")=Classes ‘tbl_df’, ‘tbl’ and 'data.frame': 5 obs. of 2 variables:
..$ INTSITE: chr "3802b" "614b" "632b" "769b" ...
..$ .rows :List of 5
.. ..$ : int 1 2 3 4 5 6 7 8 9 10 ...
.. ..$ : int 16 17 18 19 20 21 22 23 24 25 ...
.. ..$ : int 29 30 31 32 33 34 35 36 37 38 ...
.. ..$ : int 41 42 43 44 45 46 47 48 49 50 ...
.. ..$ : int 56 57 58 59 60 61 62 63 64 65 ...
..- attr(*, ".drop")= logi TRUE
unique(mercFL25miB$Species)
[1] "Bull shark " "Cobia " "Gag " "Gray snapper "
[5] "Gray triggerfish " "Greater amberjack " "Gulf flounder " "King mackerel "
[9] "Lane snapper " "Little tunny " "Red drum " "Red grouper "
[13] "Snook " "Spanish mackerel " "Vermilion snapper "
#need to convert lower case to upper case for the two dataframes to be able search properly
mercFL25miB$Species <- toupper(mercFL25miB$Species)
mercFL25miB$Species <- trimws(mercFL25miB$Species, which = c("right"), whitespace = "[ \t\r\n]")
mercFL100miB$Species <- toupper(mercFL100miB$Species)
mercFL100miB$Species <- trimws(mercFL100miB$Species, which = c("right"), whitespace = "[ \t\r\n]")
#edit COMMON SNOOK in MRIP data to read SNOOK so they match
landTtimB$common <- gsub("COMMON SNOOK", "SNOOK", landTtimB$common)
#top sites = 770, 632, 3802, 769, 614, within 25 mi
unique(mercFL25miB$INTSITE)
[1] "3802b" "614b" "632b" "769b" "770b"
unique(landTtimB$INTSITE)
[1] 770 769 632 614 3802
T614b25 <- subset(mercFL25miB,mercFL25miB$INTSITE == '614b')
T632b25 <- subset(mercFL25miB,mercFL25miB$INTSITE == '632b')
T3802b25 <- subset(mercFL25miB,mercFL25miB$INTSITE == '3802b')
T769b25 <- subset(mercFL25miB,mercFL25miB$INTSITE == '769b')
T770b25 <- subset(mercFL25miB,mercFL25miB$INTSITE == '770b')
indFL25a614 <- which((landTtimB$common %in% T614b25$Species) & (landTtimB$INTSITE == '614'))
str(indFL25a614)
int [1:5] 298 300 306 307 706
landTtimB[indFL25a614,]
unique((subset(landTtimB,landTtimB$INTSITE == '614'))$common)
[1] "PINFISH" "SCALED SARDINE" "LEATHERJACK"
[4] "SHEEPSHEAD" "ATLANTIC SPADEFISH" "STRIPED MULLET"
[7] "WHITE MULLET" "SPANISH MACKEREL" "ATLANTIC SHARPNOSE SHARK"
[10] "YELLOWFIN MENHADEN" "ATLANTIC THREAD HERRING" "GAFFTOPSAIL CATFISH"
[13] "SOUTHERN KINGFISH" "LADYFISH" "MOJARRA FAMILY"
[16] "RED DRUM" "MULLET GENUS" "FANTAIL MULLET"
unique(T614b25$Species)
[1] "BULL SHARK" "COBIA" "GAG" "GRAY SNAPPER"
[5] "GRAY TRIGGERFISH" "GREATER AMBERJACK" "GULF FLOUNDER" "KING MACKEREL"
[9] "LITTLE TUNNY" "RED DRUM" "RED GROUPER" "SNOOK"
[13] "SPANISH MACKEREL"
landTtimB$testStatus25mi <- rep('NA', nrow(landTtimB))
landTtimB$testStatus25mi[indFL25a614] <- 'Tested'
indFL25a632 <- which((landTtimB$common %in% T632b25$Species) & (landTtimB$INTSITE == '632'))
str(indFL25a632)
int [1:60] 19 20 24 25 27 28 29 30 45 52 ...
landTtimB[indFL25a632,]
unique(T632b25$Species)
[1] "BULL SHARK" "COBIA" "GAG" "GRAY SNAPPER"
[5] "GRAY TRIGGERFISH" "GREATER AMBERJACK" "GULF FLOUNDER" "LANE SNAPPER"
[9] "RED DRUM" "SNOOK" "SPANISH MACKEREL" "VERMILION SNAPPER"
landTtimB$testStatus25mi[indFL25a632] <- 'Tested'
indFL25a3802 <- which((landTtimB$common %in% T3802b25$Species) & (landTtimB$INTSITE == '3802'))
str(indFL25a3802)
int [1:65] 386 387 389 397 508 513 514 515 545 550 ...
landTtimB[indFL25a3802,]
landTtimB$testStatus25mi[indFL25a3802] <- 'Tested'
indFL25a769 <- which((landTtimB$common %in% T769b25$Species) & (landTtimB$INTSITE == '769'))
str(indFL25a769)
int [1:47] 3 6 14 16 117 122 174 212 225 226 ...
landTtimB[indFL25a769,]
landTtimB$testStatus25mi[indFL25a769] <- 'Tested'
indFL25a770 <- which((landTtimB$common %in% T770b25$Species) & (landTtimB$INTSITE == '770'))
str(indFL25a770)
int [1:15] 1 74 108 109 146 205 206 331 332 335 ...
landTtimB[indFL25a770,]
landTtimB$testStatus25mi[indFL25a770] <- 'Tested'
landTtimB$testStatus25mi[landTtimB$testStatus25mi %in% 'NA'] <- 'Not Tested'
now 100 mi, timeframe B
#top sites = 770, 632, 3802, 769, 614, within 100 mi
unique(mercFL100miB$INTSITE)
[1] "3802b" "614b" "632b" "769b" "770b"
unique(landTtimB$INTSITE)
[1] 770 769 632 614 3802
T614b100 <- subset(mercFL100miB,mercFL100miB$INTSITE == '614b')
T632b100 <- subset(mercFL100miB,mercFL100miB$INTSITE == '632b')
T3802b100 <- subset(mercFL100miB,mercFL100miB$INTSITE == '3802b')
T769b100 <- subset(mercFL100miB,mercFL100miB$INTSITE == '769b')
T770b100 <- subset(mercFL100miB,mercFL100miB$INTSITE == '770b')
indFL100a614 <- which((landTtimB$common %in% T614b100$Species) & (landTtimB$INTSITE == '614'))
str(indFL100a614)
int [1:7] 298 299 300 306 307 453 706
landTtimB[indFL100a614,]
unique((subset(landTtimB,landTtimB$INTSITE == '614'))$common)
[1] "PINFISH" "SCALED SARDINE" "LEATHERJACK"
[4] "SHEEPSHEAD" "ATLANTIC SPADEFISH" "STRIPED MULLET"
[7] "WHITE MULLET" "SPANISH MACKEREL" "ATLANTIC SHARPNOSE SHARK"
[10] "YELLOWFIN MENHADEN" "ATLANTIC THREAD HERRING" "GAFFTOPSAIL CATFISH"
[13] "SOUTHERN KINGFISH" "LADYFISH" "MOJARRA FAMILY"
[16] "RED DRUM" "MULLET GENUS" "FANTAIL MULLET"
unique(T614b100$Species)
[1] "" "ATLANTIC SHARPNOSE SHARK" "BLACKTIP SHARK"
[4] "BONNETHEAD" "BULL SHARK" "COBIA"
[7] "CREVALLE JACK" "FLORIDA POMPANO" "GAG"
[10] "GRAY SNAPPER" "GRAY TRIGGERFISH" "GREATER AMBERJACK"
[13] "GULF FLOUNDER" "KING MACKEREL" "LANE SNAPPER"
[16] "LEMON SHARK" "LITTLE TUNNY" "RED DRUM"
[19] "RED GROUPER" "SNOOK" "SOUTHERN KINGFISH"
[22] "SPANISH MACKEREL" "VERMILION SNAPPER"
landTtimB$testStatus100mi <- rep('NA', nrow(landTtimB))
landTtimB$testStatus100mi[indFL100a614] <- 'Tested'
indFL100a632 <- which((landTtimB$common %in% T632b100$Species) & (landTtimB$INTSITE == '632'))
str(indFL100a632)
int [1:66] 19 20 24 25 27 28 29 30 45 52 ...
landTtimB[indFL100a632,]
unique(T632b100$Species)
[1] "" "BLACKTIP SHARK" "BONNETHEAD" "BULL SHARK"
[5] "COBIA" "FLORIDA POMPANO" "GAG" "GRAY SNAPPER"
[9] "GRAY TRIGGERFISH" "GREATER AMBERJACK" "GULF FLOUNDER" "KING MACKEREL"
[13] "LANE SNAPPER" "LEMON SHARK" "LITTLE TUNNY" "RED DRUM"
[17] "RED GROUPER" "SNOOK" "SOUTHERN KINGFISH" "SPANISH MACKEREL"
[21] "VERMILION SNAPPER"
landTtimB$testStatus100mi[indFL100a632] <- 'Tested'
indFL100a3802 <- which((landTtimB$common %in% T3802b100$Species) & (landTtimB$INTSITE == '3802'))
str(indFL100a3802)
int [1:67] 386 387 389 397 508 513 514 515 545 550 ...
landTtimB[indFL100a3802,]
landTtimB$testStatus100mi[indFL100a3802] <- 'Tested'
indFL100a769 <- which((landTtimB$common %in% T769b100$Species) & (landTtimB$INTSITE == '769'))
str(indFL100a769)
int [1:48] 3 6 14 16 117 122 174 212 225 226 ...
landTtimB[indFL100a769,]
unique((subset(landTtimB,landTtimB$INTSITE == '769'))$common)
[1] "SAND SEATROUT" "SPANISH MACKEREL" "SCALED SARDINE"
[4] "PINFISH" "ATLANTIC THREAD HERRING" "SPOTTED SEATROUT"
[7] "WHITE GRUNT" "PIGFISH" "CREVALLE JACK"
[10] "BLACK SEA BASS" "BLUE RUNNER" "SPOTTAIL PINFISH"
[13] "GULF FLOUNDER" "GRAY SNAPPER" "BLUEFISH"
[16] "COWNOSE RAY" "GAG" "SAND PERCH"
[19] "ATLANTIC BUMPER" "ATLANTIC CROAKER" "SHEEPSHEAD"
[22] "SEATROUT GENUS" "LADYFISH" "LITTLE TUNNY"
[25] "LANE SNAPPER" "SOUTHERN KINGFISH" "SPANISH SARDINE"
[28] "NORTHERN KINGFISH"
unique(T769b100$Species)
[1] "" "BLACKTIP SHARK" "BONNETHEAD" "BULL SHARK"
[5] "COBIA" "FLORIDA POMPANO" "GAG" "GRAY SNAPPER"
[9] "GRAY TRIGGERFISH" "GREATER AMBERJACK" "GULF FLOUNDER" "KING MACKEREL"
[13] "LANE SNAPPER" "LEMON SHARK" "LITTLE TUNNY" "RED DRUM"
[17] "RED GROUPER" "SNOOK" "SOUTHERN KINGFISH" "SPANISH MACKEREL"
[21] "VERMILION SNAPPER"
landTtimB$testStatus100mi[indFL100a769] <- 'Tested'
indFL100a770 <- which((landTtimB$common %in% T770b100$Species) & (landTtimB$INTSITE == '770'))
str(indFL100a770)
int [1:16] 1 74 108 109 139 146 205 206 331 332 ...
landTtimB[indFL100a770,]
landTtimB$testStatus100mi[indFL100a770] <- 'Tested'
landTtimB$testStatus100mi[landTtimB$testStatus100mi %in% 'NA'] <- 'Not Tested'
FLbProp25mi <- landTtimB %>%
group_by(INTSITE, testStatus25mi) %>%
summarise(n_distinct(common), landings = sum(landing), biomass = sum(wgt_ab1))
FLbProp100mi <- landTtimB %>%
group_by(INTSITE, testStatus100mi) %>%
summarise(n_distinct(common), landings = sum(landing), biomass = sum(wgt_ab1))
write.csv(FLbProp25mi, "~/FUIteam/PydioData/env/data_outputs/FLbProp25mi.csv")
write.csv(FLbProp100mi, "~/FUIteam/PydioData/env/data_outputs/FLbProp100mi.csv")
Size Comparison try it out, LA time frame A
Going to try this a different way, ignoring intsite but also removing duplicate mercury results. Will only use mercury samples from within 100 mi
On to FLorida!
landTtimAwt <- landTtimA
landTtimAwt$aveWt <- (landTtimAwt$wgt_ab1*1000)/landTtimAwt$landing #in g
str(DFfl100miA)
'data.frame': 938 obs. of 37 variables:
$ Group : int 121 120 120 120 120 120 120 120 120 120 ...
$ Bottle : chr "TB04010" "ch0441" "ch0440" "TB04009" ...
$ Location.x : chr "TBM" "chm" "chm" "TBX" ...
$ Field_no : chr "TBM04110906" "chm04071102" "chm04071505" "TBX04062101?" ...
$ Spec_no : int NA NA NA NA NA NA NA NA NA NA ...
$ FIM_spec_no : chr "1" "2" "1" "1" ...
$ Other_field_no: chr "" "" "" "" ...
$ Date : int 20041104 20040712 20040709 20040621 20040616 20040616 20040615 20040613 20040613 20040613 ...
$ New_date : chr "11/4/2004" "7/12/2004" "7/9/2004" "6/21/2004" ...
$ Year : int 2004 2004 2004 2004 2004 2004 2004 2004 2004 2004 ...
$ Scfname : chr "Rachycentron canadum" "Sciaenops ocellatus" "Rachycentron canadum" "Carcharhinus leucas" ...
$ Species : chr "Cobia " "Red drum " "Cobia " "Bull shark " ...
$ NODCCODE : num 8.84e+09 8.84e+09 8.84e+09 8.71e+09 8.84e+09 ...
$ Sex : chr "F" "M" "F" "F" ...
$ SL : chr "1030" "597" "1070" "." ...
$ FL : chr "1122" "." "1170" "." ...
$ TL : chr "1279" "704" "1335" "2800" ...
$ Whole_weight : chr "." "." "20411.00" "." ...
$ Zone : chr "b" "e" "a" "" ...
$ Grid : chr "." "341" "15" "." ...
$ Lat_deg : chr "." "26" "26" "." ...
$ Lat_min : chr "." "29.374" "58.004" "." ...
$ Site : chr "big bayou" "" "" "found at gandy ramp" ...
$ Region.x : chr "MW" "" "" "MW" ...
$ Hg : num 0.8 0.22 1.5 1.9 0.26 0.23 0.35 3.1 3.5 1.9 ...
$ Lat_Deg : int 27 26 26 27 26 26 26 26 26 26 ...
$ Lat_Min : num 44.3 29.4 58 53.5 48.9 ...
$ Lon : int 82 82 82 82 82 82 82 82 82 82 ...
$ Lon_Deg : int -82 -82 -82 -82 -82 -82 -82 -82 -82 -82 ...
$ Lon_Min : num 38.38 10.14 11.17 31.98 9.01 ...
$ Region.y : chr "TB" "CH" "CH" "MW" ...
$ Location.y : chr "big bayou" "" "" "found at gandy ramp" ...
$ LatDD : num 27.7 26.5 27 27.9 26.8 ...
$ LonDD : num -82.6 -82.2 -82.2 -82.5 -82.2 ...
$ INTSITE : num 770 770 770 770 770 770 770 770 770 770 ...
$ LonDD.1 : num 136935 183158 181514 147519 185068 ...
$ LatDD.1 : num 377392 238863 291732 394346 274897 ...
unique(DFfl100miA$INTSITE)
[1] 770 769 632 615 742
subDFfl100miA <- DFfl100miA[,1:34]
str(subDFfl100miA)
'data.frame': 938 obs. of 34 variables:
$ Group : int 121 120 120 120 120 120 120 120 120 120 ...
$ Bottle : chr "TB04010" "ch0441" "ch0440" "TB04009" ...
$ Location.x : chr "TBM" "chm" "chm" "TBX" ...
$ Field_no : chr "TBM04110906" "chm04071102" "chm04071505" "TBX04062101?" ...
$ Spec_no : int NA NA NA NA NA NA NA NA NA NA ...
$ FIM_spec_no : chr "1" "2" "1" "1" ...
$ Other_field_no: chr "" "" "" "" ...
$ Date : int 20041104 20040712 20040709 20040621 20040616 20040616 20040615 20040613 20040613 20040613 ...
$ New_date : chr "11/4/2004" "7/12/2004" "7/9/2004" "6/21/2004" ...
$ Year : int 2004 2004 2004 2004 2004 2004 2004 2004 2004 2004 ...
$ Scfname : chr "Rachycentron canadum" "Sciaenops ocellatus" "Rachycentron canadum" "Carcharhinus leucas" ...
$ Species : chr "Cobia " "Red drum " "Cobia " "Bull shark " ...
$ NODCCODE : num 8.84e+09 8.84e+09 8.84e+09 8.71e+09 8.84e+09 ...
$ Sex : chr "F" "M" "F" "F" ...
$ SL : chr "1030" "597" "1070" "." ...
$ FL : chr "1122" "." "1170" "." ...
$ TL : chr "1279" "704" "1335" "2800" ...
$ Whole_weight : chr "." "." "20411.00" "." ...
$ Zone : chr "b" "e" "a" "" ...
$ Grid : chr "." "341" "15" "." ...
$ Lat_deg : chr "." "26" "26" "." ...
$ Lat_min : chr "." "29.374" "58.004" "." ...
$ Site : chr "big bayou" "" "" "found at gandy ramp" ...
$ Region.x : chr "MW" "" "" "MW" ...
$ Hg : num 0.8 0.22 1.5 1.9 0.26 0.23 0.35 3.1 3.5 1.9 ...
$ Lat_Deg : int 27 26 26 27 26 26 26 26 26 26 ...
$ Lat_Min : num 44.3 29.4 58 53.5 48.9 ...
$ Lon : int 82 82 82 82 82 82 82 82 82 82 ...
$ Lon_Deg : int -82 -82 -82 -82 -82 -82 -82 -82 -82 -82 ...
$ Lon_Min : num 38.38 10.14 11.17 31.98 9.01 ...
$ Region.y : chr "TB" "CH" "CH" "MW" ...
$ Location.y : chr "big bayou" "" "" "found at gandy ramp" ...
$ LatDD : num 27.7 26.5 27 27.9 26.8 ...
$ LonDD : num -82.6 -82.2 -82.2 -82.5 -82.2 ...
meltDFfl100miA <- subDFfl100miA %>% distinct()
str(meltDFfl100miA)
'data.frame': 280 obs. of 34 variables:
$ Group : int 121 120 120 120 120 120 120 120 120 120 ...
$ Bottle : chr "TB04010" "ch0441" "ch0440" "TB04009" ...
$ Location.x : chr "TBM" "chm" "chm" "TBX" ...
$ Field_no : chr "TBM04110906" "chm04071102" "chm04071505" "TBX04062101?" ...
$ Spec_no : int NA NA NA NA NA NA NA NA NA NA ...
$ FIM_spec_no : chr "1" "2" "1" "1" ...
$ Other_field_no: chr "" "" "" "" ...
$ Date : int 20041104 20040712 20040709 20040621 20040616 20040616 20040615 20040613 20040613 20040613 ...
$ New_date : chr "11/4/2004" "7/12/2004" "7/9/2004" "6/21/2004" ...
$ Year : int 2004 2004 2004 2004 2004 2004 2004 2004 2004 2004 ...
$ Scfname : chr "Rachycentron canadum" "Sciaenops ocellatus" "Rachycentron canadum" "Carcharhinus leucas" ...
$ Species : chr "Cobia " "Red drum " "Cobia " "Bull shark " ...
$ NODCCODE : num 8.84e+09 8.84e+09 8.84e+09 8.71e+09 8.84e+09 ...
$ Sex : chr "F" "M" "F" "F" ...
$ SL : chr "1030" "597" "1070" "." ...
$ FL : chr "1122" "." "1170" "." ...
$ TL : chr "1279" "704" "1335" "2800" ...
$ Whole_weight : chr "." "." "20411.00" "." ...
$ Zone : chr "b" "e" "a" "" ...
$ Grid : chr "." "341" "15" "." ...
$ Lat_deg : chr "." "26" "26" "." ...
$ Lat_min : chr "." "29.374" "58.004" "." ...
$ Site : chr "big bayou" "" "" "found at gandy ramp" ...
$ Region.x : chr "MW" "" "" "MW" ...
$ Hg : num 0.8 0.22 1.5 1.9 0.26 0.23 0.35 3.1 3.5 1.9 ...
$ Lat_Deg : int 27 26 26 27 26 26 26 26 26 26 ...
$ Lat_Min : num 44.3 29.4 58 53.5 48.9 ...
$ Lon : int 82 82 82 82 82 82 82 82 82 82 ...
$ Lon_Deg : int -82 -82 -82 -82 -82 -82 -82 -82 -82 -82 ...
$ Lon_Min : num 38.38 10.14 11.17 31.98 9.01 ...
$ Region.y : chr "TB" "CH" "CH" "MW" ...
$ Location.y : chr "big bayou" "" "" "found at gandy ramp" ...
$ LatDD : num 27.7 26.5 27 27.9 26.8 ...
$ LonDD : num -82.6 -82.2 -82.2 -82.5 -82.2 ...
unique(meltDFfl100miA$Species)
[1] "Cobia " "Red drum " "Bull shark "
[4] "Southern kingfish " "Bonnethead " "Blacktip shark "
[7] "" "Lemon shark " "Florida pompano "
[10] "Gulf flounder " "Red grouper " "Gag "
[13] "Snook " "Gray snapper " "Spanish mackerel "
[16] "King mackerel " "Little tunny " "Vermilion snapper "
[19] "Lane snapper " "Gray triggerfish " "Greater amberjack "
[22] "Crevalle jack " "Atlantic sharpnose shark " "Spotted seatrout "
#need to convert lower case to upper case for the two dataframes to be able search properly
meltDFfl100miA$Species <- toupper(meltDFfl100miA$Species)
meltDFfl100miA$Species <- trimws(meltDFfl100miA$Species, which = c("right"), whitespace = "[ \t\r\n]")
#now compare mrip to merc for time frame A in Louisiana
mripFLTestedA <- subset(landTtimAwt,landTtimAwt$testStatus100mi == 'Tested')
mercFLtestA <- meltDFfl100miA[meltDFfl100miA$Species %in% (unique(mripFLTestedA$common)),]
unique(mripFLTestedA$common)
[1] "SPANISH MACKEREL" "SOUTHERN KINGFISH" "SPOTTED SEATROUT"
[4] "KING MACKEREL" "FLORIDA POMPANO" "GAG"
[7] "GRAY SNAPPER" "GULF FLOUNDER" "CREVALLE JACK"
[10] "RED DRUM" "LITTLE TUNNY" "ATLANTIC SHARPNOSE SHARK"
[13] "LANE SNAPPER" "RED GROUPER" "SNOOK"
[16] "BONNETHEAD" "GRAY TRIGGERFISH"
unique(mercFLtestA$Species)
[1] "RED DRUM" "SOUTHERN KINGFISH" "BONNETHEAD"
[4] "FLORIDA POMPANO" "GULF FLOUNDER" "RED GROUPER"
[7] "GAG" "SNOOK" "GRAY SNAPPER"
[10] "SPANISH MACKEREL" "KING MACKEREL" "LITTLE TUNNY"
[13] "LANE SNAPPER" "GRAY TRIGGERFISH" "CREVALLE JACK"
[16] "ATLANTIC SHARPNOSE SHARK" "SPOTTED SEATROUT"
str(mripFLTestedA)
'data.frame': 403 obs. of 20 variables:
$ YEAR : chr "2005" "2005" "2005" "2005" ...
$ REG_RES : int 7 7 7 7 7 7 7 7 7 7 ...
$ ST_RES : int 12 12 12 12 12 12 12 12 12 12 ...
$ CNTY_RES : int 103 103 57 103 103 103 57 57 103 57 ...
$ ST : chr "12" "12" "12" "12" ...
$ CNTY : int 103 103 103 103 103 103 103 103 103 103 ...
$ INTSITE : int 615 769 770 742 770 770 770 770 770 770 ...
$ MODE_F : int 1 1 1 3 1 1 1 1 1 1 ...
$ MODE_FX : chr "3" "3" "3" "3" ...
$ ID_CODE : chr "1297920050326028" "1297920050106028" "1523120050430020" "1297920050111007" ...
$ ZIP : chr "33710" "33705" "33611" "33759" ...
$ common : chr "SPANISH MACKEREL" "SOUTHERN KINGFISH" "SPANISH MACKEREL" "SPOTTED SEATROUT" ...
$ WAVE : int 2 1 2 1 2 2 2 2 2 2 ...
$ SP_CODE : num 8.85e+09 8.84e+09 8.85e+09 8.84e+09 8.85e+09 ...
$ harv_d3 : num 0 0 0 0 0 0 0 1 1 0 ...
$ wgt_ab1 : num 3.779 1.147 0.585 0.803 2.277 ...
$ landing : num 9.04 4.02 1.01 1 6.01 ...
$ testStatus25mi : chr "Tested" "Not Tested" "Tested" "Not Tested" ...
$ testStatus100mi: chr "Tested" "Tested" "Tested" "Tested" ...
$ aveWt : num 418 285 580 800 379 ...
str(mercFLtestA)
'data.frame': 236 obs. of 34 variables:
$ Group : int 120 120 120 120 120 120 120 120 120 120 ...
$ Bottle : chr "ch0441" "ch0438" "ch0439" "ch0437" ...
$ Location.x : chr "chm" "chm" "chm" "chm" ...
$ Field_no : chr "chm04071102" "chm04060303" "chm04061702" "chm04060504" ...
$ Spec_no : int NA NA NA NA NA NA NA NA NA NA ...
$ FIM_spec_no : chr "2" "1" "1" "1" ...
$ Other_field_no: chr "" "" "" "" ...
$ Date : int 20040712 20040616 20040616 20040615 20040608 20040526 20040520 20040520 20040520 20040504 ...
$ New_date : chr "7/12/2004" "6/16/2004" "6/16/2004" "6/15/2004" ...
$ Year : int 2004 2004 2004 2004 2004 2004 2004 2004 2004 2004 ...
$ Scfname : chr "Sciaenops ocellatus" "Sciaenops ocellatus" "Menticirrhus americanus" "Sphyrna tiburo" ...
$ Species : chr "RED DRUM" "RED DRUM" "SOUTHERN KINGFISH" "BONNETHEAD" ...
$ NODCCODE : num 8.84e+09 8.84e+09 8.84e+09 8.71e+09 8.84e+09 ...
$ Sex : chr "M" "F" "F" "M" ...
$ SL : chr "597" "522" "220" "501" ...
$ FL : chr "." "." "." "." ...
$ TL : chr "704" "630" "260" "885" ...
$ Whole_weight : chr "." "." "." "." ...
$ Zone : chr "e" "a" "c" "b" ...
$ Grid : chr "341" "104" "137" "118" ...
$ Lat_deg : chr "26" "26" "26" "26" ...
$ Lat_min : chr "29.374" "48.890" "46.743" "47.601" ...
$ Site : chr "" "" "" "" ...
$ Region.x : chr "" "" "" "" ...
$ Hg : num 0.22 0.26 0.23 0.35 0.23 0.15 0.12 0.15 0.14 0.12 ...
$ Lat_Deg : int 26 26 26 26 26 26 26 26 26 26 ...
$ Lat_Min : num 29.4 48.9 46.7 47.7 47.6 ...
$ Lon : int 82 82 82 82 82 82 82 82 82 82 ...
$ Lon_Deg : int -82 -82 -82 -82 -82 -82 -82 -82 -82 -82 ...
$ Lon_Min : num 10.14 9.01 4.22 16.23 15.44 ...
$ Region.y : chr "CH" "CH" "CH" "CH" ...
$ Location.y : chr "" "" "" "" ...
$ LatDD : num 26.5 26.8 26.8 26.8 26.8 ...
$ LonDD : num -82.2 -82.2 -82.1 -82.3 -82.3 ...
mercFLtestA
AsubFLmrip <- mripFLTestedA[,c(11:12,20)]
AsubFLmrip$source <- 'MRIP'
colnames(AsubFLmrip) <- c("ZIP", "spName", "aveWt_g","source")
AsubFLmerc <- mercFLtestA[,c(12,18)]
AsubFLmerc$source <- 'MERC'
AsubFLmerc$ZIP <- 'NA'
AsubFLmerc <- AsubFLmerc[,c("ZIP", "Species", "Whole_weight","source")]
colnames(AsubFLmerc) <- c("ZIP", "spName", "aveWt_g","source")
compFLa <- rbind(AsubFLmrip,AsubFLmerc)
ggplot(compFLa, aes(x=spName, y=aveWt_g, fill = source)) +
geom_boxplot() +
xlab("Species") +
ylab("Average Weight (g)") +
theme(axis.text.x = element_text(size = 12, face = "bold", angle = 45, hjust = 1),
axis.text.y = element_text(size = 12, face = "bold"),
axis.title.x = element_text(size = 16, face = "plain"),
axis.title.y = element_text(size = 16, face = "plain")) +
theme(legend.text=element_text(size = 12, face = "plain"),
legend.title=element_text(size = 14, face = "bold"))
#not enough observations of mercury weights to run
results.FLa <-ddply(compFLa,.(spName),function(x) summarize(x,
Pvalue=t.test(aveWt_g~source,data=x,na.rm=TRUE,paired=F)$p.value,
Estimate.MERC=t.test(aveWt_g~source,data=x,na.rm=TRUE,paired=F)$estimate[1],
Estimate.MRIP=t.test(aveWt_g~source,data=x,na.rm=TRUE,paired=F)$estimate[2]))
argument is not numeric or logical: returning NAError in t.test.default(x = "117.00", y = "2860", na.rm = TRUE, paired = FALSE) :
not enough 'x' observations
Add a new chunk by clicking the Insert Chunk button on the toolbar or by pressing Ctrl+Alt+I.
When you save the notebook, an HTML file containing the code and output will be saved alongside it (click the Preview button or press Ctrl+Shift+K to preview the HTML file).
The preview shows you a rendered HTML copy of the contents of the editor. Consequently, unlike Knit, Preview does not run any R code chunks. Instead, the output of the chunk when it was last run in the editor is displayed.